Excitations of attractive 1-D bosons: Binding vs. fermionization 
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The stationary states of few bosons in a one-dimensional harmonic trap are investigated throughout the 
crossover from weak to strongly attractive interactions. For sufficient attraction, three different classes of states 
emerge: (i) iV-body bound states, (ii) bound states of smaller fragments, and (iii) gas-like states that fermionize, 
that is, map to ideal fermions in the limit of infinite attraction. The two-body correlations and momentum spec- 
tra characteristic of the three classes are discussed, and the results are illustrated using the soluble two-particle 
model. 
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I. INTRODUCTION 

In recent years, ultracold atoms have become a flexible tool 
for the simulation of fundamental quantum systems 01 Ht] . 
Their versatility derives mainly from the fact that both their 
external forces and atomic interactions can be designed to a 
great extent. One striking example is the possibility to con- 
fine the atoms' motion to lower dimensions, such as in the 
one-dimensional (ID) Bose gas yl]. Since, pictorially speak- 
ing, particles moving on a line cannot move around each other, 
they are in a sense more strongly correlated than their higher- 
dimensional counterparts. Moreover, their effective interac- 
tion strength can be tuned freely so as to enter the interest- 
ing regime of strong interactions, either via Feshbach reso- 
nances of the 3D scattering length ]4[] or through confinement- 
induced resonances of the effective ID coupling |@]. 

The case of repulsive interactions has long received con- 
siderable attention, mostly for the striking feature that, in the 
hard-core limit of infinite repulsion between the bosons, the 
system maps to an ideal Fermi gas @]. In this fermioniza- 
tion limit, the bosons become impenetrable, which has a simi- 
lar effect as Pauli's exclusion principle for identical fermions. 
The seminal exact solutions derived for special systems at ar- 
bitrary interaction strength — such as the homogeneous Bose 
gas on a ring in the thermodynamic limit (71 18[] as well as for 
definite particle numbers [9], and for the inhomogeneous gas 
in a hard-wall trap I Kill — have recovered this borderline case 
in the limit of infinite coupling. The thermodynamic nature 
of the fermionization crossover from weak to strong repulsion 
has first been explored and contrasted with the complementary 
Thomas-Fermi regime ITU [l2tl . and its microscopic mecha- 
nism has been unraveled by recent numerically exact studies 
lfl3i [Til f]~5l [l6tl . Moreover, its experimental demonstration 
has sparked renewed interest in ID bosons fl7l[T8ll . 

By contrast, the understanding of the attractive case is more 
patchy. In the homogeneous system, the ground state forms 
an iV-body bound state IU9I1 . For sufficient attraction, it be- 
comes ever more localized and is unstable in the thermody- 



namic limit. For finite systems, though, the ground state re- 
mains stable for arbitrary finite attraction, as is demonstrated 
by the exact solution via Bethe's ansatz for a ring [9] and a 
hard-wall trap ifloll . Much less is known about excited states. 
In the homogeneous system again, Monte Carlo simulations 
have indicated the existence of a highly excited gas-like state 
for ultrastrong attraction, which can be seen as the counter- 
part of the fermionized ground state for repulsive interactions 
fed. I2H1 . The evidence for this super-Tonks gas has been 
supported by a Bethe-ansatz solution l22ll . Still, an intuitive 
understanding of these states from a microscopic perspective 
and how they come about in the crossover from weak inter- 
actions is still missing. The complementary crossover for the 
low-lying excitation spectrum in turn has been investigated re- 
cently for the homogeneous system [23]; however, it does not 
include the gas-like super-Tonks gas. 

In this article, we study the entire crossover from the non- 
interacting to the strongly attractive limit for few bosons in a 
harmonic trap. This is done via the numerically exact multi- 
configuration time-dependent Hartree method introduced in 
Sec.|II] Sectionlllllpresents the general Bose-Fermi map valid 
for the gas-like super-Tonks states, and illustrates its meaning 
on the simple model of two bosons. The numerical investiga- 
tion of the stationary states in Sec. |IV] reveals three distinct 
classes for strong enough attraction: iV-body bound states 
(Sec. IIV At . states involving smaller fragments (Sec. HVBt . 
and finally gas-like states that fermionize (Sec. lIVCi i. 



II. MODEL AND COMPUTATIONAL METHOD 

Model We consider N trapped bosons described by the 
Hamiltonian 



H 



N 

£ 

i=l 



P 2 i + U( Xi ) 



i<j 



We will focus on the case of harmonic confinement, U (x) 
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x (where harmonic-oscillator units are employed through- 
out.) The effective interaction resembles a ID con- 
tact potential, but is mollified with a Gaussian S a (x) = 

-x 2 /2a 2 



I \j2/ko~ (of width a = 0.05) for numerical reasons 
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(cf. lfl4ll for details.) We concentrate on attractive forces 
g 6 (—00, 0], which can be achieved experimentally by either 
having negative scattering lengths or by reducing the trans- 
verse confinement length aj_ = y/h/Mu±_ sufficiently 

Computational method Our approach relies on the nu- 
mericallyexact multi-configuration time-dependent Hartree 
method [24l|25ll23l, a quantum-dynamics approach which has 
been applied successfully to systems of few identical bosons 
IB dl H3 IH IH as well as to Bose-Bose mixtures (H. 
Its principal idea is to solve the time-dependent Schrodinger 
equation ffi(t) = H^(t) as an initial-value problem by ex- 
panding the solution in terms of direct (or Hartree) products 
$j = <Pji ® • • • <8> <Pj„: 

$>(t) = J2Mt)®j(t). (l) 

J 

The unknown single-particle functions ipj (j = 1 , . . . , n) are 
in turn represented in a fixed basis of, in our case, harmonic- 
oscillator orbitals. The permutation symmetry of is ensured 
by the correct symmetrization of the expansion coefficients 
Aj. 

Note that, in the above expansion, not only the coefficients 
Aj but also the single-particle functions ipj are time depen- 
dent. Using the Dirac-Frenkel variational principle, one can 
derive equations of motion for both Aj, ipj B25I1 . Integrating 
this differential-equation system allows us to obtain the time 
evolution of the system via ((TJ. This has the advantage that 
the basis set {$j(t)} is variationally optimal at each time t; 
thus it can be kept relatively small. 

Although designed for time-dependent studies, it is also 
possible to apply this approach to stationary states. This is 
done via the so-called relaxation method 113111 . The key idea 
is to propagate some wave function ^(0) by the non-unitary 
e~ Hr (propagation in imaginary time.) As r — > 00, this expo- 
nentially damps out any contribution but that stemming from 
the true ground state like e -( E ™- E o)r j n practice, one relies 
on a more sophisticated scheme termed improved relaxation 
lHHHH], which is much more robust especially for excitations. 
Here (\f r |i2"|\I r ) is minimized with respect to both the coeffi- 
cients A j and the orbitals ipj. The effective eigenvalue prob- 
lems thus obtained are then solved iteratively by first solving 
for Aj with fixed orbitals and then 'optimizing' ipj by propa- 
gating them in imaginary time over a short period. That cycle 
will then be repeated. 



HI. BOSE-FERMI MAP FOR ATTRACTIVE BOSONS 

In this section, we state the general Bose-Fermi map and 
discuss its application to infinitely attractive interactions. 
(Without loss of generality, we focus on the time-independent 
formulation.) Its intuitive meaning will be illustrated on the 
special example of two harmonically trapped bosons. 



A. General map 

The Schrodinger equation of N bosons with point interac- 
tions, (E — H)^> = 0, is equivalent to a noninteracting system 
(E — Hq)^ = with boundary conditions 

2$.¥| P=0+ = s*| r=0 , (2) 

where r = Xi — Xj for fixed i ^ j. For infinitely repulsive 
interactions, g — ► 00, the constraint that 



leads to the well-known hard-core boundary conditions 

*U =!C , =0, i £ j. (3) 

Since, otherwise, fulfills the noninteracting Schrodinger 
equation, it is intelligible that it can be mapped to a nonin- 
teracting state of identical fermions, which automatically 
satisfies the hard-core boundary condition (0) by Pauli's ex- 
clusion principle [6]: 

* = A*_, A(X) = \[sgrx{x i -Xj). 

The Bose-Fermi map A serves only to restore bosonic permu- 
tation symmetry. Note that, since A 2 = 1, all local quantities 
derived from = \^\ 2 will coincide with those computed 
from the fermion state. In this sense, the case of infinite repul- 
sion is commonly referred to as fermionization limit ("Tonks 
gas"). 

By contrast, the constraint for infinite attraction 



may be satisfied by (a) d r ^\ 0+ — > —00 (provided > 
diverges more slowly than the local derivative) or (b) — > 
CP (assuming d r 1 i'\ 0+ > 0) . If case (b) applies to all i ^ j, 
then we recover the hard-core gas above. Consequently, the 
Bose-Fermi mapping ^ = A'fy- then holds for g — > —00 as 
well. In particular, the energetically lowest such state exactly 
equals thefermionized repulsive ground state - the Tonks gas. 
However, for strong attraction, this eigenstate will be highly 
excited, whereas the ground state will be strongly bound. For 
finite g > —00, this may be identified with the super-Tonks 
state in Ref. J20ll . 

A few comments are in order: 

(i) By construction, this holds for any external potential U, 
just like the standard Bose-Fermi map. In particular, the an- 
alytic solution for N bosons in a harmonic trap carries right 
over B3411 

* (X) oc e-l*! 2 / 2 \xi-xjl (4) 

l<i<j<N 

and likewise for the homogeneous system jfl . 
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(ii) By the same logic as above, this extends to binary mix- 
tures of bosons (for the repulsive case, cf. HH). Likewise, the 
generalized Bose-Fermi map for spinor bosons J36ll ought to 
apply also to the limit of infinite attraction. 

(iii) The map also holds in the presence of additional long- 
range interactions, such as in dipolar gases l37ll . 



B. Illustration 

To visualize the above argument, let us resort to the simple 
model of two bosons in a harmonic-oscillator (HO) potential, 
U(x) = t;X 2 . Here the center of mass (CM) R = J2i x i 
and the relative coordinate r — X\ — x% separate, 



H 



Pr 



1 



-NR 1 



1 



Pr 



gS{r) 



_2N 2 

One can therefore write the wave function and its energy as 

^ = 4>m®^; E = (N + i) + e, 

where <\>pj is the HO orbital with quantum number M = 
0, 1, . . . . The relative Hamiltonian may be viewed as a har- 
monic potential split into halves in the center, i.e., at the point 
of collision r = 0. There the delta function imposes the 
boundary condition (O, which amounts to a cusp for g < 0. 
The problem can be solved analytically in terms of parabolic 
cylinder functions U(a, b) 11381 13911 



where e(g) = v{g) 
dental equation 



i> e {r)=cU(-e,r), (5) 
+ | is determined through the transcen- 



u{g) G f-\0) : /» := 2 3 / 2 



r(i^) 



(6) 



r(-|) Iy 

The solution for attractive interactions is plotted in Fig.Q] 

• First imagine we start from the noninteracting ground 
state i/(0) = (Fig. [Qi). For g < 0, the Gaus- 
sian will pick up a cusp at r = 0, which becomes 
ever sharper for increasing \g\, i.e., i[i'(0 + ) — > — oo. 
In the limit g — ► — oo, when the support of ip be- 
comes much smaller than the oscillator length, this is 
described simply by the bound state of the delta poten- 



tial, ip(r) 



— 7=1 



-H/o with 



-2/g the ID scatter- 



ing length. Clearly |V>( r )| 2 9 —*°°5(r). 

An eigenstate with i/(0) = 2n ^ will also form a cusp 
at r = (Fig.QJi). In contrast to before, however, this 
peak will become ever smaller for g — > — oo, until it 
reaches down to i/j(0) — for infinite attraction. In this 
limit, the wave function's modulus equals that of the 
next lower fermionic state, ip(r) — > sgn(r)u2 n -i(r), 

where u a (r) oc H a \ -j=] e~ r2 / A . Concomitantly, 



energy 



oc 

will be 



lowered 



to lirrig^^oo v(g) 



the 

2n — 1, which matches exactly the fermionized level 

rim g _>+oo v(g) starting from P(0) = 2(n — 1). 
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Figure 1: (color online) Relative wave function ip v ( g )(r) of two 
bosons in a harmonic trap (in arbitrary units). 

(a) Ground state for attractive interactions: The Gaussian ground 
state (v(g — 0) = 0, thin-blue line) becomes peaked for g — —0.6 
(dashed-magenta). This tends to an exponential peak as g — > — oo 
(cf. g = —3.2, thick-ocher line). 

(b) Excited states for attractive interactions: The first bosonic exci- 
tation (v = 2, thin-blue line) picks up a cusp at zero for g = —1.5 
(dashed-magenta). As g — > — oo (cf. g = —16, thick-ocher), this 
cusp becomes sharper but, at the same time, is damped out more and 
more until the wave function's modulus equals the fermionic state 
ui(r) oc r e~ r 7/4 (dotted-green). 
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Figure 2: Energy spectrum {E m (g)} of N — 3 bosons in a har- 
monic trap. Different line styles correspond to states which connect 
to a trimer ( ), a dimer plus a single boson (— • — ), and a gas- 
like fermionized state ( — ), respectively. Inset: Energy E(g) for the 
lowest fermionized states of TV = 2 — 4 atoms (bottom to top). 



IV. MOLECULE FORMATION VS. FERMIONIZATION IN 
A HARMONIC TRAP 

After having presented the general mapping valid for in- 
finite attraction, let us now investigate the crossover to that 
borderline case, starting from the noninteracting states. For 
concreteness, we shall focus on the case of = 3 bosons in 
a harmonic trap. 

Figure [2] explores the evolution of the low-lying few -body 
spectrum {E m (g)} as we vary g < 0. In the absence of 
interactions, the spectrum exhibits an even spacing, which 
comes about by distributing all N particles in number states 
| AT) = \N , Nt,...) over the lowest HO levels e a = a + \. 
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Clearly, Fig.|2]reveals that, when attractive interactions g < 
are switched on, the depicted spectrum falls apart into three 
(for N = 3) qualitatively different subclasses. In anticipation 
of our analysis below, we identify the asymptotically lowest 
set of levels as TV -body bound states (or trimers), the clus- 
ter on top of these as hybrid states (dimers plus one free bo- 
son), and the highest level as gas-like state which undergoes 
fermionization. Note that, by CM separation, for each state 
there is a countable set of copies with different CM energies 
TV + |, which are shifted with respect to one another. 



A. Trimer states 

The ground state, as well as its copies with higher CM ex- 
citations, can be thought of as an iV-body bound state, which 
would be the straightforward generalization of the two-body 
bound state discussed in Sec. I11IBI In fact, this class of states 
is well known from the homogeneous system fljjjl . where an 
analytic solution is available via Bethe's ansatz fl9l l23ll 

E{g) oc -cg 2 N{N 2 -l); 9{X) oc e _ £'<* I**-"*!/", (7) 

where a = — 2/g is the ID scattering length. We will now 
argue that, for sufficient attraction, this wave function also 
holds in our case of harmonic confinement, up to a triv- 
ial CM factor 4>j^. To this end, let us proceed like in 

id 

and transform X = (x\, . . . , xn) T to Jacobian coordinates 
,r;y_i) T , here for simplicity specified for 



Y = (R,r u 
N = 3: 



Y = OX, O 
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1 
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one-body density, which represents an average over all mea- 
surements. Moreover, in contrast to the hard- wall trap i flol l. 



there is no tipping point where the width (Ax) 



c 2 ) would 



become larger again and return to its noninteracting value 
for g — ► — oo. In fact, here Ax decreases monotonically for 
g — ► — oo and eventually saturates in the CM density 



lim p(x) = \4>u(x)\ , 



(9) 



whose length scale ocm = 1/yN is suppressed by the total 
mass NM (in units of the atomic mass M = 1). To prove 
this, note that by CM separation, 



p(xx 



dxo 



• dx N \cj>x (R)\ p lci (r = {r k }). 



For increasing attraction, linig^-oo p IC \(r) = S(r), in agree- 
ment with Fig. |4] visualizes that the two-body density 
P2 {xi,X2) — (xi, Xi\p2\x\ 1 X2) is more and more concen- 
trated on the diagonal {x\ = £2}. Carrying out the trivial 
integrals proves Eq. (0, which reflects that all atoms clump 
together to point-like molecule whose position coincides with 
the CM. (Of course, the validity of our effective model re- 
quires the molecule to be still large compared to the spatial 
extension of the atoms.) 

That line of reasoning fails in the hard-wall trap, where R 
and r couple strongly due to the anharmonicity of U (x) 11401k 
in fact, for very tight binding of two particles, their common 
CM would be permitted to spread out over the whole box, thus 
compensating the stronger localization in the relative coordi- 
nate. A simple dimensional argument to see this is that the 
length scale in a hard-wall trap is simply the size of the box, 
independent of the object's mass - in contrast to the harmonic 
oscillator. 



Up to a factor, R (ri) coincide with the usual center of mass 
(two-particle relative coordinate), while T2 gives the differ- 
ence between particle #3 and the center of mass of the cluster 
(1,2). By orthogonality of O, the Hamiltonian transforms to 

H(Y) = h CM + H reU with 



2 r 



Hrel — 



E 

k=l 



1 2 , 1 2 

2 Pr > k + 2 Tk 



(8) 



If all N particles cling together to form a tightly bound state, 
then their distances will be small compared to the confine- 
ment scale, \rf-\ <C 1. In this limit, }^r\ may be safely 
neglected, so that the relative wave function asymptotically 
maps to the homogeneous form (0. Likewise, the energy 
scales as E(g) ~ Af + y — a N g 2 (a N > 0). 

A look at the one-body density profiles p(x) = (x\p\\x) 
shown in Fig. |3J a ) confirms that this state becomes more and 
more localized. In contrast to the translation-invariant homo- 
geneous case, this "soliton"-like state is localized even in the 



B. Hybrid states 

The behavior of the class of levels on top of the trimer lev- 
els (N = 3, Fig.01 is clearly more complex, which indicates 
that different CM excitations J\f are involved. On the one 
hand, the levels are significantly above those of the iV-body 
bound states; on the other hand, they do not saturate with 
increasing attraction. This suggests to identify these three- 
boson states with the formation of dimers plus one unbound 
atom. For general N, this class of hybrid states involves dif- 
ferent fragments — labeled / = 1, . . . , F — of A^-body bound 

states, where S/=i Nf — N and F = 2, . . . , iV — 1. 

Figure [3jb) shows the density profile of the lowest hybrid 
state, which has TV = 0. A pronounced peak at x — builds 
up similar to the trimer case, but with a non-Gaussian struc- 
ture indicative of relative excitations. This becomes clearer 
from a look at the two-body density in Fig. [4] For increasing 
attraction, part of the bosons clump together, as is visible from 
the "molecule" peak at {x\ = X2}. However, unlike before, 
a non-negligible part remains isolated at {a;i = — X2}. To 
understand this pattern a little better, let us revisit the Hamil- 
tonian dHJ. If we assume that only two of the three bosons 
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Figure 3: One-body density p(x) for the ground state (a), some hybrid state (b), and the lowest fermionizing state (c) of N — 3 bosons, plotted 
for different interaction strengths g. 
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Figure 4: Two-body correlation functions ^2(^1, x%) for the ground 
state, some hybrid state, and the lowest fermionizing state (from left 
to right), at coupling strengths g = —0.1, —2 and —10 (from top to 
bottom). 



bind, say |n| -C 1 (up to permutation symmetry S+), then we 
end up with two decoupled Hamiltonians, 



^.1 + ^1) 



P 2 r.2 + 



1 



9^2) 



one for a free-space molecule (ri, relating to the relative 
ground state described in Fig. QJ) and one for an effective 
particle in an excited state of an oscillator with a delta-type 
dimple at the origin (r 2 , the lowest excitation corresponding 
to Fig.QJi). This yields the relative wave function (excluding 
the trivial CM factor) 



ip{X) oc 5 4 



-\xi— x 2 \/ajj 



where the parameters are determined in analogy to Eq. ((5]). 
We have checked that this expression qualitatively reproduces 



the two-body pattern in Fig. |4]for g — > —00. That makes it 
tempting to think of the hybrid state as a hard-core "gas" of 
a dimer — clumped near the trap center — and a third, unbound 
boson. In that regime, the energy scales as E(g) ~ J\f + 
|r — 2a 5 ' + a ' wnere li m g— -00 e(g) — 5 = a € 2N — 1. 
In general, there can be different combinations of relative and 
CM excitations (TV 7^ 0) which give nearly the same energy 
- this explains the splittings of all but the lowest hybrid levels 
in Fig. |2] 



C. Fermionizing states 

Let us now focus on the highest level in Fig. [2] which is 
the energetically lowest gas-like state. Its energy does not 
diverge quadratically with g — > —00, but rather saturates. 
By the fermionization map above, its limit is simply the en- 
ergy of N free fermions, liirig^-oo E(g) = Y, a =n e & = 
(see Fig. |2j inset), and likewise for higher excitations. Ev- 
idently, this requires a huge energy for the connecting level 
E(g = 0) > E(—oo) if N > 1. In fact, the difference be- 
tween the two equals that between the noninteracting ground 
state Neo and its fermionization limit ^^=0 £ a, which can 
be written down explicitly in a harmonic trap, 

N 2 N N(N - 1) 
\E(g -> -co) - £(0)| = ^" - f = [ 2 j - 

This can be thought of as increasing (lowering) the energy 
by Ae = 1 for each pair (i < j). Therefore, the "super- 
Tonks" level connects to the noninteracting level £'(0) = 
ig_ + n(n-i) = N ^ N _ 1 ^ as may be verified for N < 4 in 

Fig- inset). 

Accordingly, the corresponding many-body state is ex- 
pected to evolve to @. A glimpse of this can be got- 
ten from the density profile shown in Fig. [3jc), where the 

density wiggles characteristic of the Tonks gas emerge, 

P = yZa=o l^«l 2 ( a similar observation has been stated in 
Ref. Il2lll ). In the repulsive case, this has the familiar interpre- 
tation that the N bosons localize on N more or less "discrete" 
spots due to a trade-off between mutual isolation and external 
confinement CUfMl. 
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Some more insight into this crossover is given by the 
two-body correlation function p?,(x\,x?) (Fig- 0- As 
expected from the two-atom toy model (Sec. IIIIBI i. the 
diagonal {x± — X2} "damps out" more and more. The 
fact that it persists even for couplings as large as g = —10 
underscores the notion of the (finite-<?) super- Tonks gas being 
more strongly correlated than its repulsive counterpart l2(ill . 
This in turn relates to the picture that, due to a positive ID 
scattering length a = —2/g, a small region is excluded from 
the scattering zone, so that the hard core effectively extends 
to a nonzero volume. This also offers an explanation for 
another phenomenon: As g — ► —00, we see that the typical 
fermionized checkerboard pattern forms in the two-body 
density lfl5l l34ll . This signifies that, upon measuring a first 
boson at, say, x\ ps 0, the remaining N — 1 = 2 bosons are 
pinpointed to discrete positions X2 ~ ±1.5. However, here 
the peaks are much more pronounced than in the Tonks gas, 
which may be accounted for by a "thicker" hard core between 
the atoms. 

We have so far looked into local observables, where the 
fermionization limit imprints truly fermionic properties on the 
bosons. By contrast, nonlocal features as evidenced, e.g., in 
the experimentally relevant momentum distribution 



p(k) =2ir{k\ Pl \k) = 




reflect nontrivial differences from the ideal Fermi gas. Fig- 
ure |5jc) shows the evolution of p(k) for the gas-like state in 
juxtaposition with the trimer and dimer states [Figs. |3a-b)]. 
For the ground state, where all bosons simply form an ever 
tighter TV-body molecule as g — ► — 00, pi(x, x') oc S(x — x') 
loses all long-range order, i.e., p\{x, x') — for \x — x'\ > 0. 
By complementarity, its momentum spectrum trivially ap- 
proaches a flat shape. Things are more complicated for the 
hybrid state in Fig.|5lb): Since only two bosons bind and, as a 
whole, form a hard-core composite with the remaining atom, 
some long-range order is preserved, so that the central peak 
p(Q) persists even for large values of \g\. Note that, as in the 
repulsive case, the hard-core short-range correlations enforce 
an algebraic decay for high momenta, p(k) ~ c/k A 14 111 . Fi- 
nally, the gas-like state exhibits the most interesting behavior 
(Fig. |5p). The initially box-like distribution p(k) = 2wp(k) 
(harmonic trap at g — 0) forms a strong k = peak highly 
reminiscent of the Tonks gas lfl3l, [l5ll . Also note the slow for- 
mation of the characteristic k~ 4 tails. This complements our 



picture of the crossover from zero to strongly attractive inter- 
actions. 



V. SUMMARY 

In this work, we have brought together the subjects of at- 
tractive, one-dimensional Bose gases — which currently are of 
great interest and experimentally relevant — and the binding 
properties of few-body systems. We have studied the sta- 
tionary states of one-dimensional bosons in a harmonic trap 
throughout the crossover from weakly to strongly attractive 
interactions. 

Three different classes of states have emerged for strong 
enough attraction: (i) The ground state and its center-of-mass 
excitations become A -body bound states, for which any two 
particles pair up to a tightly bound molecule. Its binding 
length a shrinks to zero with increasing attraction, and thus 
the relative motion becomes independent of the trap geome- 
try, (ii) By contrast, certain highly excited states fermionize, 
i.e., they map to an ideal Fermi state for infinite attraction. 
Both the typical fermionic density profile with N maxima and, 
more generally, the characteristic checkerboard pattern in the 
two-body density have been evidenced, signifying localiza- 
tion of the individual atoms. Also, the formation of a hard- 
core momentum distribution has been witnessed, with a zero- 
momentum peak and an algebraic decay for large momenta, 
(iii) Between these two extremes, there is a rich class of hy- 
brid states featuring mixed molecule and hard-core boundary 
conditions. For the special case of N = 3 atoms, this class 
consists of a dimer plus a single boson, with a hard-core sep- 
aration between the two. 

Even though we have focused mostly on few atoms (A = 
3) in a harmonic trap, these results reflect the microscopic 
mechanism for arbitrary atom numbers and external poten- 
tials. 
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